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We relate the frequency of the Josephson plasma resonance in layered superconductors with the 
frequency dependent superconducting response. We demonstrate that the sharp resonance can 
persist even when the global superconducting coherence is broken provided the resonance frequency 
is larger than the frequency of interlayer phase slips. In this situation the plasma frequency is 
determined by the average Josephson energy, which can be calculated using the high temperature 
expansion. We also find the temperature dependence of the average Josephson energy from the 
Monte Carlo simulations and determine the applicability region of the high temperature expansion. 
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Interlayer Josephson coupling between superconduct- 
ing layers has crucial influence on thermal and pin- 
ning properties of layered superconductors in the mixed 
state. A new, powerful tool to probe the Josephson cou- 
pling in a wide range of temperatures and fields has 
been introduced recently. A sharp microwave absorp- 
tion resonance has been discovered in the vortex state of 
BiiSr?,CaCuiO% 0-^, which was shown to be caused 
by electric field across the layers || and was attributed 
to the earlier predicted plasmon resonance [Q-Q. It was 
found that the resonance can be tuned by temperature 
and magnetic field. The resonance field corresponding 
to the fixed microwave frequency of 30-90 GHz increases 
with temperature at low temperatures having a maxi- 
mum at the irreversibility temperature and afterwards 
smoothly decreases with temperature. The resonance 
frequency at a given temperature and field can serve as 
a measure of the effective Josephson coupling. Unex- 
pectedly, the plasmon resonance has been observed up to 
high temperatures ^75 K where the Josephson cou- 
pling should be strongly suppressed by thermal fluctua- 
tion of pancake vortices. It was theoretically predicted 
that pancake fluctuations should actually destroy the 
global superconducting coherence above the decoupling 
line [ fl2| , [l3| . Matsuda et al || suggested that the decou- 
pling line, if it exists, should lie above the temperature- 
field line at which the sharp plasmon resonance at a given 
microwave frequency is observed. However, their line of 
resonance, Bq(T), clearly lies above the decoupling line 
observed by I — V measurements J§| and falls into the 
region where linear resistivity in the c-direction was ob- 
served |8| [TT| . In this Letter we resolve this controversy. 
We demonstrate that the sharp plasmon resonance can 
actually persist even above the decoupling transition pro- 
vided the resonance frequency is higher than the typical 
frequency of the interplane phase slips. In this situa- 
tion the plasma frequency is determined by the average 
Josephson energy E e J^ = Ej (cos <p n n+i) (or, equiva- 



lently, the average Josephson current) as was suggested 
in Refs. |J|]. Here (f> nn+ i = 4> n+ i - <j> n - ^A z is the 
gauge invariant phase difference between layers and Ej 
is the bare Josephson energy. The decoupling transition 
manifests itself by a disappearance of the linear supercon- 
ducting response along the c-axis at zero frequency while 
the average Josephson energy should not change much at 
this point and it actually stays finite at any temperature. 

The superconducting response, Q^, determines the re- 
sponse of the electric current to the external alternating 
vector-potential, EG] 
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In superconducting state Q^^o > 0. For the lay- 
ered superconductor with the Josephson energy Ej — 
—Ej J dr^2 n cos 4>nn+ii one can obtain a general ther- 
modynamic expression for the z-component of Q u us- 
ing the fluctuation-dissipation theorem (FDT)(see, e.g., 
Rcf . llffl ) . The external vector-potential, A u , creates the 
energy perturbation, £ A = J dr J2 n sin 0„„+iA^ 
with jj = ij^-Ej being the Josephson current. 
The average current along the z-direction j z — 
' ' ' ' • 1 in linear with respect to A^ 
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approximation is given by 
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Here (•■•)() notates the thermodynamic average with 
undisturbed energy and notates averaging which 

takes into account the perturbation £a- The second term 
describes the average Josephson current, which emerges 
due to the phase perturbations caused by the external 
field. Using the FDT relation |l7j this term can be con- 
nected with the correlation function of the Josephson cur- 
rents and we obtain the following expressions for the real 
{Q'u ) an d imaginary (Q") parts of the superconducting 
response along the ^-direction 
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cos(ui) — (sin 0oi (0, 0) sin cj) nn+1 (r, t)) 
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Substituting the material relation (Q) into the Maxwell 
equations one can see that the plasma resonance exists 
at the frequency u> p i connected with Q' u by relation 
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provided 3> QZ- Here e is the dielectric con- 
stant. Due to the frequency dependence of QZ the 
last expression is actually an equation for wj rather 
than its definition. As one can see from Eq.(0) Q' u is 
given by the sum of two terms, the first is determined 
by the average Josephson current, jj (cos (f> nn +i), and 
the second is determined by the correlation function of 
the Josephson currents. In the decoupled liquid phase 
QL=o — 0, which means that at u> = the second 
term in square brackets of Eq.(||) exactly compensates 
the first one. The frequency dependence is totally de- 
termined by the second correlation term. This term 
is strongly suppressed when the frequency exceeds the 
typical "phase slip" frequency oj ps , which is determined 
by the typical decay time of the correlation function 
(sin0„ n+ i(O, 0) sin0„ n+ i(O, t)). In the high frequency 



regime, to ^> uj ps , we have Q' u w ^sjj (cos< 



L ), i.e., 



the plasma resonance probes the average Josephson cur- 
rent. In the mixed state (cos0 n „+i) is suppressed by 
misalignment of pancake vortices induced by pinning and 
thermal fluctuations. In Ref. [Eij this quantity has been 
connected with the density correlation functions assum- 
ing Gaussian distribution for 4> nn +i- 

In the decoupled liquid phase (cos nra +i) can be calcu- 
lated using the high temperature expansion with respect 
to the Josephson energy which gives 
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is the phase correlation function for a single 2D layer. In 
the vortex liquid state the decay length of Si(r) is given 
by the average intervortex spacing and Eq. (^J) can be 
rewritten as 



(COS0„„ + 1) = -gjTg/' 
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where / = ^ J dr ^|/Si(r)| 2 \ is the dimensionless func- 
tion of order unity with weak temperature dependence. 
We find from simulations that the the high-T expansion 
is quantitatively accurate if (cos<f> n n+i) < 0.5. This cor- 
responds to the field and temperature range TB > QqEj. 
Combining Eqs. (||) and (^) we obtain 
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Taking into account the weak temperature dependence 
of /, this expression gives an explanation for the exper- 
imentally observed dependence uj^ cx 1/TB Note 
that such scaling indicates almost decoupled layers rather 
than a persistence of coupling at high temperatures. 

The magnitude and shape of the function f(T) de- 
pend upon the physical realization of the 2D liquid state 
in the layers. In pin free superconductors this function 
is universal and depends only on the dimensionless pa- 
rameter T(47rA a b) 2 /s ( E ) §. Below we calculate this uni- 
versal function using Monte Carlo simulations. In real 
samples random pinning increases disorder in vortex ar- 
rangements even in the liquid state, which leads to extra 
suppression of f(T) and further smoothes its tempera- 
ture dependence. On the other hand, the pinning effect 
becomes progressively weaker at higher temperatures. 

Substituting (§) into (|) we obtain the high-T expan- 
sion for QZ 
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with S(r,t) = (exp [i (<p nn+1 (r,t) - nn+ i(0, 0))]) being 
the dynamic phase correlation function. One can im- 
mediately see that, indeed, in this regime QZ=q = 0, 
indicating decoupled layers. 

The result (JsJ) can be generalized to the case when the 
magnetic field is tilted with respect to the c-axis. In pres- 
ence of the y-component of the magnetic field, B y , and in 
the lowest order with respect to Ej, the phase difference 
acquires an extra contribution 4> nn +i — > 4>nn+i + i§^ByX, 
which gives an extra factor exp(i^-B y x) in the RHS of 
Eq. (JTJ). Approximating S\(r) by a Gaussian function 
Si(r) « exp(— zfig ) we obtain 
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The Josephson coupling in tilted fields in presence of 
coexistent Abrikosov and Josephson lattices was stud- 
ied in Ref. [^5). As follows from Eq. ( pi]) the ty pical 
afe-component of the field is given by yj B z ^>q/tts 2 and 
coincides with the one found in Ref. M| at high fields 
By > $ /7s 2 . For B z = 100 G this typical field is ~ 
1.7 T which corresponds to a tilt angle ^0.3°. In con- 
trast to the lattice state, the average Josephson current 
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does not experience oscillations caused by commensura- 
bility of Abrikosov and Josephson lattices but, instead, 
decreases monotonically with B y . 

One can expect that in the liquid state the phase slips 
are strongly assisted by thermal motion of pancake vor- 
tices. In real samples this motion is strongly hindered by 
pinning as follows from the thermally activated behav- 
ior of the in-plane component of resistivity p a b [ ^6|pd] 
in a wide temperature range. An important experimen- 
tal fact is that the temperature dependencies of p a b and 
p c are characterized by the same value of the activa- 
tion energy |T^,[l8|,|l9| . This indicates that the interlayer 
phase slips are, indeed, induced by motion of pancake 
vortices. These phase slips are very well described by 
a simple model which assumes that the dissipation pro- 
cesses are caused by a small concentration n of mobile 
pancakes with a thermally activated diffusion constant 
D [fl8|| . Within this model the dynamic phase correla- 
tion function decays in time and space as 
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where r cut is the in-plane phase coherence length. Sub- 
stituting the last equation into Eq. (110) we obtain 
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where the typical phase slip frequency, oj ps , is given by 
iVps = 2nnD In (nr 2 ut ). Because the same mobile pan- 
cakes determine the dissipation along the a&-plane, this 
frequency can be connected with p a b Pqj 
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orw ps w 2-f4-10 8 [l/s]-T[K]-p ab [/iO-cm]. Taking T = 45 
K and p a b — lCP 3 /£57-cm corresponding to the resonance 



field B w 0.1T |3| we obtain uj ps w 1 2 ■ 10 7 [l/s], 
which is indeed much smaller than the microwave fre- 
quency u = 2ir ■ 4.5 • 10 10 [1/s] used in Ref . j| . As follows 
from Eqs. (|,[l|) Q" « (uj ps /uo)Q' u at w » uj ps , i.e. the 
plasma resonance is not suppressed by the dissipation 
in this regime. Eq. ([jj]) holds only in the "phase slip" 
regime, i.e., until the quasiparticle contribution to the 
c-axis transport can be neglected. 

The analytical calculation of the correlation function 
(0) in the liquid state is not possible even for a pin-free su- 
perconductor. To obtain quantitative information about 
the Josephson coupling in the liquid state and to check 
the high-T expansion we perform Monte Carlo simula- 
tions of the frustrated three dimensional XY model which 
was used to investigate phase transitions in the vortex 
state pl],E2| . The energy of the model is given by 



E = J2 V « (<M n + S a ) - 0(n) - Oa(n)) 



Here <fr(n) are the phases of the order parameter located 
at the sites of a mesh 1 < n x ,n y < N±,l < n z < N z , with 
periodic boundary conditions in all directions. The vec- 
tor potential a Q (n) (a = x, y, z ) is taken in the Landau 
gauge a y (n x ) — 2itbn X) where the dimensionless mag- 
netic field b determines the fraction of the grid sites oc- 
cupied by vortices. We chose N± = 72, N z = 4, 10, 
and b = 1/36 which gives 144 vortex lines. The in- 
plane phase-phase interaction function, V x {6) — V y (9), 
determines the in-plane phase stiffness. The choice of 
the low angle asymptotics, V x , y (0) f» — 1 + 2 /2, fixes 

the energy and temperature scale as s$q/ ^vr(47rA) 2 ^. 
An unrealistic feature of the model is that vortices are 
subject to the periodic pinning potential created by the 
numerical grid. This effect can be reduced significantly 
by optimizing the shape of V x (9). We chose V x (9) — 
-3r/4 - (1 - r) cos(6>) - (r/4) cos(26>), and adjusted the 
coefficient r to minimize the energy barrier for the vor- 
tex jump to the neighbor site, which gives r = 0.376. 
The coupling between planes is given by the standard 
Josephson interaction, V z {9) — — jtCos((9), with the am- 
plitude determined by the anisotropy parameter T. The 
coupling between the two dimensional lattices in the lay- 
ers is determined by a scaled magnetic field h = B / B cr , 
where B cr — (f>o/(7s) 2 is the crossover field |l2| and 7 is 
the anisotropy of the London penetration depth (in di- 
mensionless units h = Tb). The parameter T should not 
be considered as the anisotropy of a real layered super- 
conductor but rather as a parameter which allows us to 
vary the scaled magnetic field h. We apply the standard 
Metropolis algorithm to the Hamiltonian (|l^). The evo- 
lution of the phase diagram with changing anisotropy will 
be reported elsewhere. In this Letter we will focus on the 
calculation of the quantity (cos 1) an d verification of 
the high-T expansion for this parameter. 
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FIG. 1. Temperature dependencies of the superconducting 
response at uj — (filled symbols) and parameter (cos0 nn +i) 
(open symbols) for two values of the anisotropy parameter 
T = 100 and 225 

Fig. [l] shows the numerically obtained temperature 
dependencies of the low frequency superconducting re- 
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sponse Q' = Q' u=0 and the parameter (cos</> nn +i) for 
two values of the anisotropy (r = 100 and 225). As 
one can see from this plot the temperature depen- 
dencies of these parameters are almost identical above 
the level (cos0 lm+ i) = 0.6. Below this level Q'(T) 
starts to drop fast, indicating the decoupling transi- 
tion, and (cos0„„ + i) still decreases smoothly. For high 
anisotropies (cos</> rm +i) has a small drop at the decou- 
pling transition. One can see also that for T = 100, 
Q'(T) shows a strong size effect for n z < 10, which in- 
dicates that this parameter is not simply determined by 
the interaction of two neighbor layers but rather char- 
acterizes the superconducting coherence along the whole 
sample thickness. 

To verify the high-T expansion we calculate the phase 
correlation function (Q) at different temperatures and in- 
tegrate it to determine the universal function f(T) in 
Eq. (||). The last equation takes the following form in 
dimcnsionlcss units 



(cos </»„„+!) = f(T)/(2TTb). 



(15) 



We found that the temperature dependence of f(T) is, 
indeed, rather weak and can be approximated by the for- 
mula /(f) = 0.862 - 0.976 f for 0.1 < f < 0.6. For 
parameters corresponding to BiSSCO, this interval cov- 
ers the temperature range 24 K < T < 66 K. In Fig. ^ we 
plot the temperature dependencies of the directly calcu- 
lated r (cos^nn+i) at different values of T together with 
it's high-T expression ([TJj). As one can see from this plot, 
all curves merge with the same universal curve at high 
enough temperatures which is very well approximated 
by the high-T expansion. Systems with lower anisotropy 
merge with this curve at higher temperatures. 
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FIG. 2. Temperature dependencies of the parameter 
F (cos 4>nn+i) for different V. High-T expansion curve for this 
parameter given by Eq. (JT^) is also plotted. 

We can use these results for quantitative analysis of 
the plasma resonance data at high temperatures. Using 
the relation jj = 2c$o/ [s(477rA a b) 2 ] one can extract the 
anisotropy ratio 7 from Eq. (j9[) . Taking the resonance 
field B a = 0.065 T at T = 53.9 K and u/2ir = 45 GHz 



from Ref. |3|] and using the material parameters e = 20 
and X ab (T) = 1800l/^l - (T/T c ) 2 we calculate from 
the simulation results (cos0 nn +i) ~ 0.11 and 7 w 330. 
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